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ABSTRACT 

We report the first results from the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, 
which utilizes the pure parallel orbits of the Hubble Space Telescope to do deep imaging along a large 
number of random sightlines. To date, our analysis includes 26 widely separated fields observed by 
the Wide Field Camera 3, which amounts to 122.8 arcmin 2 in total area. We have found three bright 
Yogs-dropouts, which are candidate galaxies at z > 7.4. One of these objects shows an indication 
of peculiar variability and its nature is uncertain. The other two objects are among the brightest 
candidate galaxies at these redshifts known to date (L > 2L*). Such very luminous objects could 
be the progenitors of the high-mass Lyman break galaxis (LBGs) observed at lower redshifts (up to 
z ~ 5). While our sample is still limited in size, it is much less subject to the uncertainty caused by 
"cosmic variance" than other samples because it is derived using fields along many random sightlines. 
We find that the existence of the brightest candidate at z w 7.4 is not well explained by the current 
luminosity function (LF) estimates at z fts 8. However, its inferred surface density could be explained 
by the prediction from the LFs at z « 7 if it belongs to the high-redshift tail of the galaxy population 
at z w 7. 

Subject headings: cosmology: observations — galaxies: luminosity function, mass function 



1. INTRODUCTION 

The Wide Field Camera 3 (WFC3) on board the 
Hubble Space Telescope (HST) has begun breaking new 
grounds in our exploration of the early universe. Us- 
ing the ultra-deep observations taken by the "HUDF09" 
program (PI: Illingworth), samples of candidate galax- 
ies at z 7 have been greatly enlarged (Oesch et al. 
2010; Bunker et al. 2010; McLure et al. 2010; Yan et 
al. 2010a), and candidates have been found at z w 8 
(Bouwens et al. 2010b; McLure et al. 2010; Yan et al. 
2010a) and possibly out to z ~ 10 (Yan et al. 2010a; 
Bouwens et al. 2010c). The studies of the stellar pop- 
ulations of galaxies at z > 7 are now underway (e.g., 
Bouwens et al. 2010a; Labbe et al. 2010a, 2010b; Finkel- 
stein et al. 2010). The pencil-beam HUDF09 data are 
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augmented by the shallower but 8x wider observations 
(Windhorst et al. 2010) of the WFC3 Early Release Sci- 
ence Program 2 (ERS2; PI: O'Connell), which provide 
us a unique opportunity to study the bright end of the 
luminosity function (LF) at z > 7 because of the much 
increased sample size of bright candidates (Wilkins et al. 
2010a, 2010b; Bouwens et al. 2010d; Yan et al. 2010b). 
Even with the ERS2 data, however, we are only barely 
able to probe the bright end at L ~ L*. It is important 
to extend our investigation to the brighter regime, as it 
could be where the progenitors of the high stellar-mass 
(> 10 10 M Q ) galaxies at z 6-5 would locate (Yan et al. 
2006). Such studies will need to survey larger areas, and 
preferably along different sightlines in order to reduce 
the bias caused by the underlying large scale structures 
(a.k.a. "cosmic variance"; see, e.g., Robertson 2010). 

Our Hubble Infrared Pure Parallel Imaging Extra- 
galactic Survey, or "HIPPIES" , was designed for this 
purpose. This program utilizes the unique "pure paral- 
lel" observing mode of HST, which enables simultaneous 
observations using instruments that are not doing the 
primary observations. As the pointings of the primary 
observations are drawn from a wide variety of programs 
of different objectives, the associated parallel observa- 
tions are in effect observing completely random fields, 
and hence are ideally suited for minimizing the impact 
of cosmic variance. In HST Cycle 17, there were two 
pure parallel programs using WFC3 as the imaging in- 
strument, one led by our team (PID 11702, PI: H. Yan) 
and the other led by M. Trenti et al. (PID 11700). Our 
program has been extended to the coming Cycle 18 (PID 
12286). HIPPIES will incorporate the data from all these 
programs for its science objectives, one of which being 
addressing the very bright-end of galaxy LF at z « 7 
and beyond. Here we report our initial results based 
on the data taken by the end of Cycle 17. Throughout 
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this Letter, we use the following cosmological parame- 
ters: 9,m = 0.27, n A = 0.73 and H = 71 km s" 1 Mpc" 1 . 
The quoted magnitudes are all in the AB system. 

2. DATA DESCRIPTION 

In Cycle 17, pure parallels were only allowed when the 
prime instrument was either the Cosmic Origins Spectro- 
graph (COS) or the Space Telescope Imaging Spectro- 
graph (STIS). Programs 11700 and 11702 only requested 
the high galactic latitude opportunities at \b\ > 20°, and 
only requested long-duration opportunities (> 3 orbits) 
so that the observations could be sufficiently deep. The 
observations and data from these two programs are de- 
scribed below. 

2.1. Observations 

These two programs used both the infrared (IR) and 
the UV-optical (UVIS) channels of WFC3. The IR im- 
ages were taken in the same three bands in both pro- 
grams, namely, F098M, F125W, and F160W (hereafter 
^098 1 ^105; and #125, respectively). In the UVIS channel, 
program 11700 used F606W (V(506)i while program 11702 
used F600LP (LPq 00 ). The left panel of Figure 1 shows 
the system response curves in these passbands. 

Our current work includes 26 WFC3 pure parallel fields 
acquired by the end of Cycle 17 (2010 August). The to- 
tal exposure time in these fields ranges from 6 to 34 ks, 
with the median of 10.8 ks (<~ 5 orbits). The distribu- 
tion of exposure time in the four bands roughly follows 
the ratio of UVIS:F098M:F125W:F160W= 1:2:1:1, 
but varies significantly among the fields. The CCD ob- 
servations in UVIS always have at least two CR-SPLIT 
exposures that enable the rejection of cosmic rays. Most 
of the IR observations also have at least two exposures 
each band in each field. As the IR array was always 
non-destructively read-out multiple times during the ex- 
posure (at least five times for these observations), cosmic 
ray rejection generally is not a problem, even with single 
exposures. 

2.2. Data Reduction 

Our data reduction starts from the outputs from the 
HST On-the-Fly Reprocessing (OTFR) pipeline, which 
are photometrically calibrated images in units of count 
rate (i.e., the "*_fit.fits" images fetched from the HST 
data archive). 

The vast majority of our parallel observations are non- 
dithered because the primary COS or STIS observations 
arc not dithered. This makes it difficult to reject im- 
age defects using the usual outlier clipping method in 
the stacking process. Most of the defects are flagged in 
the "data quality" (DQ) extension of the processed im- 
age, which records the pixels that are diagnosed by the 
processing software to be problematic for various rea- 
sons (e.g., bad pixels, problems in telemetry, unstable 
response, etc.). In order to reduce their chance of con- 
taminating our dropout sample, such image defects in the 
IR were filled by interpolation of nearby pixels. Using the 
archival ERS2 data, we also identified several thousands 
of additional outlier pixels that deviate from linear re- 
sponse at 5<7 level, and these pixels were also fixed in 
the same way. The defected pixels were not treated for 
the UVIS images in the current work, because these im- 
ages are only used for the "veto" purpose in the dropout 



selection, and hence such image defects in UVIS do not 
likely cause any contamination. 

The MultiDrizzle software (Koekemoer et al. 2002) 
was then used to correct for the geometric distortions 
and to stack the images in each band. For each field, 
the World Coordinate System (WCS) of the first im- 
age in V098 was always chosen as the reference to align 
all images. The final pixel scale that we adopted is 
0'.'09 pixel" 1 . The pixel noise of the drizzle-combined 
images is correlated at small scale because of the sub- 
pixel sampling. We followed the procedure utilized in 
the GOODS program (Dickinson ct al. 2004) to calcu- 
late the correlation amplitudes in our mosaics, and then 
to derive the so-called "RMS maps" that register the ab- 
solute root-mean-square noise in each pixel. 

2.3. Photometry 

Matched-aperture photometry was carried out by run- 
ning the SExtractor program of Bcrtin & Arnouts (1996) 
in dual-image mode. The ,/125-band mosaics were used 
as the detection images. We adopt the SExtractor 
MAG_AUTO magnitudes, which were measured using 
the default (Kron factor, minimum radius) of (2.5, 3.5). 
A 2-pixel, 5x5 Gaussian filter was used to convolve the 
Ji25 _ band images for detection; the detection threshold 
was set to 1.5(7, and a minimum of four connected pix- 
els above the threshold was required. The RMS maps 
were used to measure the noise in the pixels within 
the MAG_AUTO apertures. We only include sources at 
S/N> 5 in the Ji 25 -band (within the MAG.AUTO aper- 
ture) for further analysis. 

The depth of our mosaics varies significantly because 
the exposure time and the background condition vary sig- 
nificantly among the HIPPIES fields. For reference, the 
median 5er depths measured within a O'.'2-radius aperture 
are 27.36, 27.08, 27.16, 27.05, and 26.68 mag in U 606 , 
LPqoo, logs, J125, and Hi 60 , respectively. The source 
count at J125 < 25.2 mag is complete in every field, and 
is complete at J125 < 26.0 mag for 20 out of the 26 fields. 

3. CANDIDATE GALAXIES AT Z > 7.4 

In this section, we describe the selection of candidate 
galaxies at z > 7.4 as lo98-dropouts using HIPPIES data. 

3.1. Y"o98- dropout selection 

Because Y"o98 does not overlap with J125, it is possible 
to adopt a large color decrement such that the selection 
is less prone to the confusion with the 4000 A break of red 
galaxies at low redshifts and less affected by photometric 
errors. Following Yan et al. (2010b), our main color cri- 
teria are Y"o98 — ^125 > 1-3 mag and non-detections (i.e., 
S/N< 2) in the "veto" U 6 06 and/or LP 600 images. Based 
on our simulation using a large number of galaxy tem- 
plates from the models of Bruzual & Chariot (2003) and 
the line-of-sight H I absorption recipe of Madau (1995), 
the adopted lo98 — ^125 color decrement is appropriate 
for selecting galaxies at z > 7.4. As our UVIS images 
arc not very deep, and we do not have data at longer 
wavelength (such as Spitzer IRAC data; see, e.g. Yan et 
al. 2010a, 2010b), we further require that a legitimate 
candidate should have J125 — -Hi60 < 0.3 mag in order 
to further reduce the chance of contamination from red 
galaxies at low redshifts. The right panel of Figure 1 
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Fig. 1. — (left) System response curves of the WFC3 passbands 
used in the current pure parallel observations. Both programs used 
the same filters in the IR. In the UVIS, program 11700 used Vgoe 
while program 11702 used LP§qq. The vertical dotted line indicates 
the location of Lyct line redshifted to z = 7.4. (right) Color-color 
diagram demonstrating the selection of Yogs-dropouts. The black 
dots are the field objects. The shaded area indicates the selection 
area, where the two objects in our final sample are shown in solid 
red squares and the peculiar Yggg-dropout that has an indication of 
variability is shown in open red square. The connected circles to left 
show the color track of a typical young galaxy at high redshifts us- 
ing the model from Bruzual & Chariot (2003 (BC03); r = 0.3 Gyr 
and age of 100 Myr), and attenuated by the intervening H I ab- 
sorption along the sightline using the formalism of Madau (1995). 
Different colors are used to code different redshift ranges along 
these tracks: blue is for z < 6.4, magenta is for 6.4 < z < 7.3, and 
black is for z > 7.4. The selection area is well separated from the 
regions occupied by the potential sources of contaminations. The 
yellow stars show the colors of Galactic brown dwarfs (e.g. Leggett 
et al. 2002). The green symbols show the colors of a typical red 
galaxy at z f» 1-3 simulated using a BC03 model (r = 50 Myr 
and age of 2.0 Gyr). The cyan symbols represent the colors of the 
same red galaxy but with a strong [OIII] emission line of restframe 
equivalent width of lOOA. The black and the green arrows represent 
the reddening vectors in these passbands appropriate for z = 7.4 
and z = 2, respectively, using the reddening law of Calzetti (2001) 
and E(B - V) = 0.5 mag. 
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shows the selection of Yo98-dropouts on the J125 
versus Y098 — ^125 color-color diagram. 

The selected candidates were visually inspected to ex- 
clude objects that were formally reported to have S/N< 2 
in the UVIS image but still visible to eyes. In this step, 
we also rejected objects that were caused by artifacts 
and image defects that are not identified in the data re- 
duction process (see Section 2.2). A last concern about 
the contamination is that the interpolation of the identi- 
fied defected pixels might leave some residuals that could 
mimic the colors of Yo98-dropouts by chance. To deal 
with this problem, we created masks of the interpolated 
pixels (see Section 2.2) for the individual images, and 
ran MultiDrizzle to project them to the same reference 
frame of the science images. The locations of the candi- 
dates were examined on these projected masks, and we 
only kept the objects that do not have any defects in the 
central 6x6 pixel area. 

3.2. Sample 

Our selection above resulted in a total of three lo98 _ 
dropouts over 26 HIPPIES fields that amount to 122.8 
arcmin 2 in area. The colors of these three Yo98 _ dropouts 
are superposed in the right panel of Figure 1 (solid red 
squares). Their image are shown in Figures 2 and 3, and 
Table 1 lists their photometric information. 

One of these three objects, YsDrop03 (see Figure 3), 
shows an indication of peculiar variability, which makes 
its candidacy as a high-redshift galaxy less certain. We 
performed photometry in each individual image, and ob- 





Fig. 2. — Image stamps of the two Yogg-dropouts m our sample. 
Images are 5'' 4 X 5'.'4 in size. The red circles are 0'.'5 in radius. 
North is up and East is to left. 
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Fig. 3. — (top) Same as Fig. 2, but for the variable object 
(YsDrop03). The field was observed by both programs 11700 and 
11702, and hence has UVIS data in both V606 an d Weoo- (bot- 
tom) Light curves of HIPPIES YsDrop03. The J125 and H 16 o 
data points are matched-aperture MAG-AUTO magnitudes ob- 
tained using the master J125 stack as the detection image, while 
the arrows are the 2cr upper limits measured within an r = $.'2 
aperture. The errors are ltr error within the MAG-AUTO aper- 
ture. 



tained its light curves as shown in Figure 3. The variabil- 
ity is the most obvious in Hiqq, where it is the brightest 
in the first image, declines by ~ 1.0 mag within 2 days, 
and remains nearly constant in the following ^1.2 days. 
Similar trend is also seen in J125, albeit less significantly. 
However, such a variability does not seem to present in 
Y098 or the two UVIS bands, where this object always 
remains undetected although the observations in these 
bands were executed both before and after the first im- 
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Fig. 4. — Comparison of the surface density of Y098 -dropouts at 
the very bright end inferred from our sample (red solid pentagons) 
and the predictions from various LFs at z ss 7 and 8 (curves). The 
inset zooms in the region at the very bright end. The 2 7 LFs 
are taken from Bouwens et al. (2008; Bo08), Ouchi et al. (2009; 
Ou09), McLure et al. (2010; M10), Ocsch et al. (2010; OelO), 
Bouwens et al. (2010d; BolOd), and Yan et al. (2010b; Y10b), 
while the 2 fs 8 LFs are taken from BolOd, M10 and Yan ct al. 
(2010a; YlOa). The 2 si 7 LFs are integrated over 6.4 < z < 7.7, 
while the z 8 ones are integrated over 7.4 < z < 9.4. The surface 
density of galaxies at z 7 inferred from the 2850-dropouts in the 
HUDF (black open squares; Yan et al. 2010a) and the ERS2 field 
(black solid squares; Yan et al. 2010b) are also plotted. 

ages were taken in H±eo and J125. Therefore, this object 
still remains a legitimate Yo98 _ dropout after this scrutiny. 
If it is not at z > 7.4, the only viable explanation for its 
color and its variability is that it is a flaring brown dwarf 
star within our Galaxy (e.g., Schmidt et al. 2007). How- 
ever, this interpretation is less favored, because this ob- 
ject has full width at half maximum (FWHM) of ~ 0'.'3 
in both J125 and H\qo stacks and hence seems to be re- 
solved. Another possibility is that this object is an active 
galactic nucleus (AGN; see, e.g., Cohen et al. 2006) at 
z > 7.4, or a transient at a similar redshift. Without 
further data, it is difficult to determine its nature. In 
any case, it does not fit in the conventional picture of 
an LBG, and therefore we do not include it in our final 
sample. 

The remaining two are good candidate galaxies at 
z > 7.4. None of them show time variability. Both are 
resolved (FWHM > 0'/3-0'.'4) in J 125 and H 160 bands 
and hence are not likely brown dwarfs. Finally, they are 
not likely red galaxies at z « 1-3 whose 4000 A break 
could mimic Lyman break, as their J125 — i?i60 colors 
are much bluer than those of such interloppers. This 
is particularly true for YsDrop02, which has a very blue 
color of J125— ti\m ~ —0.2 mag. A few similar cases have 
already been observed in candidate galaxies at z « 7 and 
beyond (e.g., Oesch et al. 2010; Bouwens et al. 2010a; 
Yan et al. 2010a, 2010b). 

4. DISCUSSION 



The two Yogs-dropouts in our final sample are the 
brightest candidate galaxies at z > 7 discovered by 
WFC3 to date. Even assuming that their redshifts are 
at the lower end of the selection window (z « 7.4), we 
obtained their absolute magnitudes of M — —21.6 ±0.16 
and —21.2 ± 0.15 mag for YsDropOl and YsDrop02, re- 
spectively, which are at least 2x brighter than any cur- 
rent estimates of the M* value at z R3 7 or 8. Earlier 
surveys over contiguous fields have resulted in a number 
of z-band dropouts (candidate galaxies at z ~ 7) of com- 
parable brightness (e.g., Ouchi et al. 2009; Hickey et al. 
2010; Castellano et al. 2010) or even brighter (Capak et 
al. 2009), and our current work indicates that such high 
luminosity galaxies could exist at z > 7.4. 

Using the formalism of Madau et al. (1998), the 
rest-frame UV luminosities of our candidates corre- 
spond to star formation rates (SFRs) of 23.7+ 3 'i and 

16.4^2'i Mq yr -1 . As an order-of-magnitude estimate, 
if they keep forming stars at the same rates to z ~ 5, 
they could accumulate stellar masses of 0.7-1.1 x 10 10 Af Q 
in the next ~ 473 million years. Therefore, our Yogs- 
dropouts could be the progenitors of the high-mass LBGs 
observed at z ~ 5 that have stellar masses to the order of 
- 10 10 M Q (see e.g., Yan et al. 2005; Stark et al. 2007). 

The existing LF estimates of LBGs at z ~ 7 and be- 
yond are largely based on observations at L < L*. For 
this reason, it is interesting to see whether the num- 
ber density of such very bright candidate galaxies at 
z > 7.4 is consistent with the expectations from such 
LFs. This comparison is shown in Figure 4, where the 
result from this study is shown as the red pentagons. 
As the formal redshift selection window of Y"o98-dropout 
spans 7.4 < z < 9.4, we compare to the predictions from 
the LFs at both z 8 and 7. Our data point is appar- 
ently higher than the predictions from the z ks 8 LFs, in- 
cluding the one of BolOd that predicts the highest surface 
density. If we single out the brightest object YsDropOl, 
the discrepancy is more obvious. The BolOd z w 8 LF 
predicts a cumulative surface density of ~ 2.7 x 10~ 3 per 
arcmin 2 to < 25.6 mag over 7.4 < z < 9.4, which is a 
factor of three lower than the density inferred from our 
YsDropOl (8.1 x 10~ 3 per arcmin -2 ). On the other hand, 
the LF of Yan et al. (2010b; M* = -20.33, a = -1.80 
and $ = 5.52 x 10~ 4 Mpc~ 3 ) predicts the highest bright- 
end surface density at z « 7. This LF, like others that 
are mostly based on the WFC3 zsso-dropout results, is 
applicable over the redshift range of6.4<z<7.7. It 
predicts a cumulative surface density of ~ 1.5 x 10~ 2 per 
arcmin 2 to < 25.6 mag, which seems to be able to explain 
the existence of YsDropOl if it is at z < 7.7. 

5. SUMMARY 

In this Letter, we present the initial results from the 
HIPPIES program. To date, our analysis includes 26 
widely separated _ffSTWFC3 pure parallel fields at |6| > 
20° obtained in Cycle-17, which amount to 122.8 arcmin 2 
in total area. Using these data, we search for candidate 
galaxies at z > 7.4 as Yo98-dropouts. One candidate 
that we found shows an indication of peculiar variability 
in J125 and -H160 but remains undetected in the bluer 
bands, and its nature is unclear. Excluding this object, 
our current sample consists of two very bright candidates 
at L > 2L* , which, based on their SFR estimates, could 
be linked to the progenitors of high-mass LBGs observed 
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TABLE 1 
Properties of Yo98-DROPOUTS a 



ID 


R.A. and Decl. (J2000) b 


^606 


LPqoo 


^098 


Jl25 


#160 


Y-J 


J-H 


HIPPIES- YsDropOl 
HIPPIES- YsDrop02 
HIPPIES- YsDrop03 c 


16:31:35.24 +37:36:14.03 
14:36:50.56 +50:43:33.69 
07:50:51.41 +29:16:17.54 


> 28.1 

> 28.4 

> 28.4 


> 28.3 


27.33+0.71 

> 28.4 

> 29.0 


25.52+0.16 
25.86+0.15 
26.62+0.17 


25.33+0.21 
26.07+0.27 
26.46+0.20 


1.81 

> 2.5 

> 2.4 


0.19 
-0.21 
0.16 



a ' All magnitude limits arc 2<r limits measured within an r — 0. 2 aperture. 

The quoted coordinates are based on the astrometric solutions provided by the OTFR pipeline for the first Yo98 image in each field. 
c ' This object shows indication of variability and thus its nature is unclear. 



at lower redshifts (up to z ~ 5) . While its size is still very 
limited, our sample is constructed from a large number of 
random fields and thus the impact of cosmic variance is 
minimal. The surface density inferred from our brightest 
candidate is not well explained by the existing LF esti- 
mates at z ~ 8, but could be explained by the prediction 
of LFs at z « 7 if it belongs to the high-redshift tail of the 
z « 7 galaxy population. Our new HIPPIES data soon 
to be taken in the coming HST cycle (~ 42 pointings) 
will enable us to search for different yios-dropouts in a 
similar way, and the comparison of the two results will 
shed light to the evolution of the star formation activities 
in the early universe. 
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